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Reduced equations that describe low-frequency plasma dynamics play an important role in our 
understanding of plasma behavior over long time scales. One of the oldest paradigms for reduced 
plasma dynamics involves the ponderomotive Hamiltonian formulation of the oscillation-center dy- 
namics of charged particles (over slow space-time scales) in a weakly-nonuniform background plasma 
perturbed by an electromagnetic field with fast space-time scales. These reduced plasma equations 
are derived here by Lie-transform and variational methods for the case of a weakly-magnetized 
background plasma. In particular, both methods are used to derive explicit expressions for the pon- 
deromotive polarization and magnetization, which appear in the oscillation-center Vlasov-Maxwell 
equations. 
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I. INTRODUCTION 

Our understanding of complex plasma phenomena can be significantly improved by the asymptotic elimination (or 
decoupling) of fast space-time scales from slow space-time scales of interest. Over the past 30 years, Lie-transform 
perturbation methods [H, 0, Q have played an important role in the process of dynamical reduction in plasma physics. 
Standard examples of reduced plasma dynamics include the elimination of fast space-time scales associated with 
the gyromotion of a charged particle in a strong magnetic field (i.e., the guiding- center reduction 0, [H, [f| and the 
gyrocenter reduction 0, Hf) or the elimination of the fast eikonal space-time scales associated with charged-particle 
motion perturbed by a hi gh- frequency, short- wavelength electromagnetic wave (the oscillation-center reduction ft, 

HE HI 111 IE HE IE HE 111)- 

Note that, while the process of dynamical reduction is implicitly associated with the elimination of fast space- 
time scales from a set of dynamical equations, the process may or may not explicitly reduce the number of degrees 
of freedom. When it does (e.g., guiding-center reduction), a fast orbital angle (used to describe particle motion) 
becomes asymptotically ignorable and its conjugate action becomes an adiabatic invariant. Whether or not the 
number of degrees of freedom is reduced, however, the numerical integration of these dynamically-reduced equations 
of motion over long space-time scales of interest, which generally allows for more realistic plasma geometries to be 
considered, presents an important practical application in plasma physics. 

Variational formulations of reduced plasma dynamics exist for single-particle Hamiltonian dynamics as well as 
Vlasov-Maxwell (kinetic) and fluid self-consistent theories. First, the variational formulation for reduced single- 
particle dynamics relies on the dynamical reduction of the phase-space Lagrangian by Lie-transform methods |3j. 
The reduced single-particle dynamics is expressed in terms of a reduced Hamiltonian and a reduced Poisson-bracket 
structure (derived from the symplectic part of the reduced phase-space Lagrangian). Conservation laws are associated 
with the invariance of the reduced Hamiltonian on fast degrees of freedom (e.g., ignorable or fast angles) and the 
corresponding invariants (e.g., exact or adiabatic actions) appear explicitly in the symplectic part of the reduced 
phase-space Lagrangian. Second, the variational formulation of the reduced Vlasov-Maxwell (kinetic) equations, on 
the other hand, is intimately connected to the variational formulation of reduced single-particle dynamics through 
a constrained variational principle in extended phase space Conservation laws for the reduced Vlasov-Maxwell 
equations are derived by applying the Noether method on the reduced Lagrangian density. Third, the variational 
formulation for reduced fluid dynamics relies on the existence of the reduced fluid Lagrangian density, from which 
exact conservation laws are derive once again by the Noether method [2(| ■ 

A. Berkeley School of Plasma Physics 

In all of these developments, the Berkeley School of Plasma Physics, led by Allan N. Kaufman in collaboration 
with his graduate students and postdocs, has played a fundamental role. Allan's pioneering work in plasma theory 
over 50 years (with over 100 papers), which is reviewed in his memoirs (to appear in the KaufmanFest Proceedings), 
can be divided into three separate periods. During the first period (from 1957 to 1969), Allan's work focussed on 
plasma dynamics in a magnetic field (in collaborations with S. Chandrasckhar and K. M. Watson) and the statistical 
physics of an imperfect gas. During the second period (from 1970 to 1987), Allan's work focussed on nonlinear plasma 
waves, the Hamiltonian formulation of plasma physics (which included the Lie-transform perturbation analysis of 
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oscillation-center and guiding-center dynamics), the construction of dissipative Poisson brackets (by combining an 
entropy-conserving Poisson bracket and an energy-conserving dissipative bracket), wave & particle chaos, and wave 
kinetic theory. During the third period (from 1987 to the present), Allan's work focussed on linear wave conversion 
in plasmas and fluids, as well as the non-eikonal formulation of wave-action conservation laws in plasma physics. 

The main characteristics of the Berkeley School of plasma physics, which were developed mainly during the second 
and third periods of Allan's work, involve the development and applications of Hamiltonian and Lagrangian methods 



in plasma physics. In fact, with a few exceptions (e.g., R. L. Dewar [21|, P. J. Morrison [22j, and J. A. Krommes 
|23|), the use of these methods in plasma physics is intimately associated with Allan's name. It is worth mentioning 
that the Berkeley School has also greatly benefited from the influx of many postdocs (from Princeton's PPPL and 
elsewhere), visitors, and collaborators (see Allan's memoirs for details). 

The development and applications of Hamiltonian and Lagrangian methods in linear wave conversion (during the 
third period) are reviewed by Tracy and Brizard in a separate paper in the KaufmanFest Proceedings. Here, I present 
the derivation of the oscillation-center dynamics of charged particles (developed during the second period) in a weakly- 
magnetized background plasma perturbed by high-frequency electromagnetic field fluctuations. The emphasis will be 
placed on applications of Lie-transform perturbation methods (with John Cary 0] and Robert Littlejohn Q, two of 
Allan's graduate students, acting as major participants in their developments) and applications of variational (action) 
principles to the development of the ponderomotive Hamiltonian oscillation-center theory (with Allan's students John 



Cary and Bruce Boghosian [Vi 
(To| , and Philippe Similon [3, El 



as well as a series of Allan's postdocs including Shayne Johnston [9j , Celso Grebogi 
16j acting as major participants in their developments). The purpose of the paper 
is not to present a complete survey of the historical development of oscillation-center Hamiltonian dynamics in plasma 
physics (which is briefly discussed in Allan's memoirs). Instead I wish to present a new (and hopefully interesting) 
application of Lie-transform perturbation and variational methods in plasma physics. 

In the present paper, the ponderomotive (eikonal-averaged) polarization and magnetization effects due to a high- 
frequency electromagnetic wave propagating in a weakly-magnetized plasma are investigated by two complementary 
methods. The ponderomotive polarization and magnetization are derived either from an oscillation-center vari- 
ational principle involving the weak background electromagnetic fields (Eo,Bo) as variational fields, or from the 
(Lie-transform) push-forward relation between the particle fluid moments and the oscillation-center fluid moments. 
An important consequence of the variational formulation is that an exact energy-momentum conservation law can be 
derived by the Noether method and is expressed in terms of linear and ponderomotive polarization and magnetization. 
We note that while some of the fundamental issues raised by this work are resolved within a covariant relativistic treat- 
ment, we present a non-relativistic treatment here (except when discussing the moving-magnetic-dipole controbution 
to polarization in Sec. IIII D[) and postpone the relativistic treatment for future work. 



B. Organization 

The remainder of the paper is organized as follows. In Sec. [H] the variational formulation of the exact Vlasov- 
Maxwell equations is reviewed [l9j . The formulation is given in terms of extended phase-space coordinates, in which 
the energy-time canonical-coordinate pair is added to the six-dimensional regular phase-space coordinates. From the 
Vlasov-Maxwell action functional, we also derive conservation laws by using the Noether method. In Sec. IIIII the 
variational formulation of a generic set of reduced Vlasov-Maxwell equations is presented [24[ • Here, the terms "exact" 
and "reduced" are used to emphasize the fact that, while the exact Vlasov-Maxwell equations exhibit the full range 
of space-time scales associated with exact particle Hamiltonian dynamics, the reduced Vlasov-Maxwell equations 
exhibit only the long space-time scales associated with reduced Hamiltonian dynamics. The explicit derivations of 
general expressions for the polarization and magnetization by push-forward (Lie-transform) method and by variational 
method are presented. The energy-momentum conservation laws for the self-consistent reduced Vlasov-Maxwell 
equations are also derived by the Noether method. In Sec. II V[ the Lie-transform perturbation analysis of oscillation- 
center Hamiltonian dynamics is presented in terms of non-canonical coordinates, so that both the Poisson-bracket 
(symplectic) structure and the Hamiltonian are perturbed by the electromagnetic wave fields. Next, the linear and 
nonlinear polarization and magnetization are derived by variational and push-forward methods. Lastly, a summary 
of this work is presented in Sec. [V] 
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II. EXACT VLASOV-MAXWELL EQUATIONS 
A. Extended Hamiltonian Dynamics 

The most general setting for Hamiltonian perturbation theory [24| is the eight-dimensional extended phase space 
with coordinates z a = (t, x; w, p), where w and p denote the kinetic energy-momentum coordinates and (x, t) denote 
the space-time location of a charged particle (of mass m and charge e). Here, the Hamiltonian description of charged- 
particle dynamics in an electromagnetic field F Ml/ = d^A v — d v A^, represented by the four-potential = (— $, A) 
[using the space-time metric g^ v — diag(— 1, 1, 1, 1)], is expressed in extended phase space in terms of (i) the extended 
Hamiltonian 

H(z) = J- |p| 2 - w = 0, (1) 
Ira 

and (ii) the extended phase-space Lagrangian (summation over repeated indices is implied) 

r = (p + -A) -dx - (w + e$) dt = + dx" = T a dz a , (2) 

where x M = (ct, x) and p^ = (— w/c, p). The constraint TL = in Eq. (TTJ) implies that the physical motion takes place 
on the surface w = \p\ 2 /2m. 

The extended equations of motion are obtained from the variational principle 

= 5 - Hdr^j , (3) 

where r is the Hamiltonian orbit parameter in extended phase space. The variational principle ([3]) yields the Eulcr- 
Lagrange equations 

dz b dH 

where the Lagrange matrix uj (with components uj a b = d a Tb — dbT a ) is associated with the differential two-form 



uj = dr = dp^ A drr/ 1 + — dx 1 " A dx" = - u ab dz a A dz b . (5) 



The extended Poisson bracket { , } is obtained from the extended phase-space Lagrangian © by inverting the 
Lagrange matrix (with components uj a b) to obtain the Poisson matrix J = u)^ 1 . From the Poisson-matrix components 
J ab = {z a , z b }, we obtain the extended Poisson bracket 



1 ' 1 \dx» dp„ d Pfl dx*) c ^ dp„ d Vv dz a u dz b7 



(6) 



defined in terms of two arbitrary functions F and G. Hamilton's equations in extended phase space are expressed as 

£ - H) = J-W (7) 



which includes 

dx^ _ on 
dt dp^ 



= V* = (c, v = , (8) 



and dp^/dt = (e/c) F^v", when dt/dr = 1 is substituted. 

Next, the extended Vlasov equation, which describes the time evolution of the particle distribution on extended 
phase space, is expressed in terms of the extended Hamiltonian (fl]) and extended Poisson bracket © as 

dT dz a dTU) , 

dr dr dz a 1 J w 
where, in order to satisfy the physical constraint ([1]), the extended Vlasov distribution is defined as 

T(z) = c5(w-\p\ 2 /2m)f(x,p,t), (10) 
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and /(x, p, t) denotes the time-dependent Vlasov distribution on regular phase space. By integrating the extended 
Vlasov equation ^ over the energy coordinate w (and using dr = dt), we obtain the regular Vlasov equation 

= df = d£ dx d£ dp df_ 

dt ~ dt dt ' dx dt ' dp' [ ' 

Lastly, the extended Vlasov equation © is coupled to Maxwell's equations 

V-E = 4vrp, (12) 

VxB--— - — J, 13 

c at c 

where E = V$ - cr 1 dA/dt and B = V X A satisfy the constraints V • B = and V X E + c _1 <9 t B = 0, and the 
charge-current densities 

are defined in terms of the extended Vlasov distribution (p~0|) . with d A p = c~ 1 dwd 3 p. 



B. Vlasov-Maxwell variational principle 

The Vlasov-Maxwell equations §§§ and (|T5|) - (jT5]) can be derived from the Vlasov-Maxwell action functional [l9| 

A[T,A, t ] = - E / d & zT(z) H(z) + J ^ F : F = J d 4 x C(x). (15) 
The variation of the Lagrangian density 

C(x) = ±- F(x) : F(x) - E / d *P H*,P) ( 16 ) 

yields (after rearranging terms) 

SC , SA^l 9 ^ + £ - E W ^ + (17) 

To obtain this expression, we used the constrained (Eulerian) variation for T, defined as 

ST = AT - Sz a d a T = - 6z a d a T, (18) 

where the Lagrangian variation AT = (by definition since the Vlasov distribution is constant along a Lagrangian 
orbit in phase space). The virtual displacement in extended phase-space is 

Sz a = - ({S, z a } + -J {x^z a }) , (19) 

where S is the canonical generating scalar field for this phase-space displacement. Hence, by substituting the dis- 
placement (fl"9)) into Eq. |[T8|) . we obtain the Eulerian variation 

e dT 

5T = {S, T} + -5A^—. (20) 
Lastly, the last term in Eq. (fTT)) involves the Noether four-vector 

A" -i*v- + £/ IS- (2D 

which does not contribute to the variational principle 

5Cd i x = (22) 
because this term appears as a space-time divergence in Eq. (|17[) . 
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1. Vlasov-Maxwell Equations. 

We obtain the extended Vlasov equation ([9]) by simply requiring stationarity in Eq. (|22|) with respect to S. Sta- 
tionarity with respect to 5A V , on the other hand, yields the Maxwell equations 

£*P._-£J/^|?^ (23, 

Decomposition in terms of components of F yields Eqs. (|12 | - (fT3|) . Note that the electromagnetic field tensor also 
satisfies the Maxwell constraint equations d a V ^ + d^V va + d v V all = 0. 

2. Energy-Momentum Conservation Law. 

An important advantage of a variational formulation is that one can use Noether's Theorem [25} to derive ex- 
plicit conservation laws. The conservation laws for the Vlasov-Maxwell equations © and (f2"5| are derived from the 
Noether equation SC = d u K" . By substituting explicit expressions for the variations (S,SA^,SC), we can obtain the 
conservation laws of energy-momentum, angular momentum, and wave action. The latter conservation law plays a 
crucial role in the analysis of linear wave conversion in plasmas (see paper by Tracy and Brizard in the KaufmanFest 
Proceedings). 

We now present the derivation of the energy-momentum conservation law associated with invariance of the La- 
grangian density ([To]) with respect to space-time translation x^ — ► x^ + Sx^. Under this translation, the variations 
(S, SA^.SC) are 

S = (p M + eA u /c) Sx^ ) 
5A^ ee F^ 8x v — (A u 6x v ) }, (24) 
SC ee ~d u (5x^C) J 

where S generates the virtual displacement Sx^ = {x 11 , S} as required. After substituting the variations (f2~4")l in 
the Noether equation pTjl , and using the Maxwell equations (|23|) , we obtain the energy- momentum conservation law 
<9 M T MI/ = 0, where the stress-energy tensor 



-jp.il ^ _}_ 
4ir 



i^"(F:F) - (F-Fr 



is the sum of the field contribution (involving the tensor F) and the particle Vlasov contribution (involving the Vlasov 
distribution J 7 ). 



III. REDUCED PLASMA DYNAMICS 



The processes of dynamical reduction in single-particle plasma dynamics and plasma kinetic theory are associated 
with the extended near-identity phase-space transformation T e : z — > z = 7^z, and its inverse T" 1 :z^z = 7^ _1 z. 
These transformations are expressed as asymptotic expansions in powers of a small dimensionless ordering parameter 
e representing either a space-time scale ordering or the strength of the perturbation (e.g., wave amplitude) as 0] 



e G\ + e 2 (G§ + \ d • dGf) 
e Gf - e 2 (G§ - | Gx ■ dGf) 



(26) 



Here, the nth-order vector field G„ • d = G b n db is used to eliminate the fast space-time scales at order e n . Note that 
in standard applications of Lie-transform perturbation theory, the time coordinate is unchanged (i.e., t ee t) so that 
G l n ee at all orders. 

A. Reduced Push-forward and Pull-back Operators 



The near-identity transformation (|26|) induces a transformation of the extended Vlasov equation (|9J) as follows. First, 
the push- forward operator T~ x associated with the inverse phase-space transformation Xr 1 , defined in Eq. (|26|) . leads 
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to the scalar-invariance relation 

F(z) = m = T(T-^) = T-^(z) (27) 

between the particle Vlasov distribution T on the extended particle phase space and the reduced Vlasov distribution 
T on the extended reduced phase space. 

Next, we construct the reduced extended Vlasov operator d e / 'dr, defined in terms of the push- forward operator T7 1 
(and its inverse, the pull-back operator T e ) acting on an arbitrary function G as 



dr ~ e 



= {O, H},. 



In the last expression, d e /dr is expressed in terms of the extended reduced Hamiltonian TL and the extended reduced 
Poisson bracket 

[F, g} £ = T7 1 ({T e F, T,g}) . (28) 

The reduced Poisson bracket ([28]) can also be constructed from the reduced phase-space Lagrangian u) e = dT e = 
d (T^T + dS) = T-^dr) = I^uj, so that jf = (u^ 1 )"'' = ^h, where we used the fact that d 2 S = (i.e., 
V X = in three dimensions) and the exterior derivative d commutes with the operators T e and T^ 1 - 

Using the reduced extended Vlasov operator d e /dr induced by the near- identity phase-space transformation (|26l) . 
the extended reduced Hamilton's equations are expressed as 

{2°, n} e , (29) 



dr 

where the extended reduced Hamiltonian 

H = T^H = H-w (30) 

is defined as the push-forward of the extended Hamiltonian H = H — w, with H = T~ 1 H — dS/dt. For many 
practical applications of reduced plasma dynamics, however, the new energy-momentum coordinates p^ = (— w/c,p M ) 
are canonical coordinates and, therefore, the reduced Poisson bracket { , } e = { , } c is the canonical bracket 

{f, d} c = (dP/dx^) (dd/dp^) ~ (dP/dp^) (dd/a®*). 

Lastly, the near-identity phase-space transformation (|26l) also introduces the reduced-displacement vector 

p e = T7 x x - x, (31) 

defined as the difference between the push- forward T^x of the particle position x and the reduced position x. The 
reduced displacement ([31]) is expressed as 

p e = -eG? - c^GJ-iGi-dGf) +••■ (32) 

in terms of the generating vector fields (Gi, G2, • • •) associated with the near-identity phase-space transformation (|26[) . 
This reduced displacement plays in important role in the expression of polarization and magnetization effects in 
reduced plasma dynamics (see Sec. MI E)]) . 

B. Reduced Vlasov-Maxwell Equations 

The extended reduced Vlasov equation for the extended reduced Vlasov distribution F is expressed as 

= \T, H \ = = (33) 

I ' Jc dr dr dz a v ; 

Here, the reduced Vlasov distribution is defined as the push-forward of the particle Vlasov distribution 



F(z) = c5[w-H(x,p,t)]F(%p,t), 



(34) 
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with the reduced extended Hamiltonian (|30|) satisfying the physical constraint H = 0. When integrated over the 
reduced energy coordinate w, Eq. (33) becomes the regular reduced Vlasov equation 

d e F _ dF rf e x dF d£ dF 

dt ~ dt dt ' 5x dt ' dp' 1 ' 

where we used dt/dr = 1 to eliminate r. 

The dynamical reduction associated with the phase-space transformation (|26[) introduces polarization and magne- 
tization effects into the Maxwell equations, which transforms the microscopic Maxwell's equations (fT^ - tflTJl) into the 
macroscopic (reduced) Maxwell's equations [U HI] 

V-D(x) = ^irp(x), (36) 

i an 4ir - 

VxHx - -^-(z) = -J4 (37) 
cot c 

Here, the terms "microscopic" and "macroscopic" are used to emphasize the fact that, while the sources (p, J) of the 
microscopic Maxwell equations fT2")) - (fl"3")) are expressed in terms of moments of the particle Vlasov distribution 
the sources (p, J) of the macroscopic Maxwell equations (|36|) - (|37p are expressed in terms of moments of the reduced 
Vlasov distribution T . In Eqs. (|5oT) - ([3T|) , the microscopic electric and magnetic fields E and B are replaced by the 
macroscopic fields [27j 



D(x) = E(ar) + 47rP(a;) 
H(a;) = B(jc) - 47rM(a;) 



(38) 



where P and M are the polarization and magnetization associated with the dynamical reduction (|26p . In the present 
paper, these polarization and magnetization are derived either by variational method (based on the existence of a 
reduced Lagrangian density) or by push-forward (Lie-transform) method (which establishes the connection between 
fluid moments in particle phase space and fluid moments in reduced phase space). 

By comparing the microscopic and macroscopic (reduced) Maxwell's equations, we note that the dynamical reduc- 
tion associated with the phase-space transformation ([26]) has introduced the following expressions for the charge and 
current densities: 

p(x) = p(x) - V-P(x), (39) 

3<x) = 3(x) + ^-(x) + cVxM(x), (40) 
ot 

where p po i = — V-P denotes the polarization density, J po i = dP/dt denotes the polarization current, and J mag = 
cV X M denotes the magnetization current. Note that the reduced charge-current densities (p, J) satisfy the reduced 
charge conservation law 

f + v-5 = t + |(v-r) + v. J -v.(f + cvx, 

= f + V.J S 0. (11) 

The macroscopic (reduced) Maxwell's equations (f36"|) -([37 j) can also be written in terms of the microscopic fields (E, B) 

as 

V-E = 4tt(p - V-P), (42) 

1 <9E 4tt /- dP \ 
VXB--- = - J+-+CVXM , (43) 



where the polarization charge density and polarization-magnetization currents appear explicitly. It is this microscopic 
form that is most often useful for practical applications of reduced plasma dynamics [20(. 
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C. Reduced Vlasov-Maxwell Variational Principle 



We now show that the reduced Vlasov-Maxwell equations (|33[) and (|42|) - (I43|) can be derived from the reduced 
variational principle J d A x 8C = 0, where the reduced Lagrangian density is 

C(x) ee F(x) : F(x) - E / d 4 p T(x,p) H(x,p; A,F). (44) 

Note that, as a result of the dynamical reduction of the Vlasov equation, the reduced Hamiltonian appearing in 
Eq. (I44p is expressed in terms of canonical energy-momentum coordinates as 



H(x,p;A,F) ee ^—\$--A(x)\' + e$(x)-w) + * e (p - - A(x); F(x) J , (45) 

where the reduced ponderomotive potential ^ depends explicitly on the field tensor F liv {x). From these field depen- 
dences, we define the reduced four-current density 

> = (C p,j) = -c E / - E e / ( 46 ) 

where 

c?i dp^ \ ' m VI ~ c 7 dp 
The reduced antisymmetric polarization-magnetization tensor [17], on the other hand, is defined as 



d e x» an ( i _ e A . , 

c, - (p - - A) + — . (47) 



where the polarization and mag netization K° l = P l and K ij = e ijfe M fe are defined as 



Reduced Vlasov-Maxwell Variational Principle. 



By following the same steps outlined in Sec. Ill Bl we begin our variational formulation of the reduced Vlasov-Maxwell 
equations with an expression for the variation of the Lagrangian density 



5C = - E / di p£ ' {P, n} 



6A V 

47T 



-in 



r 



where the reduced Noether four-vector 

A^ EE E / ""/'^ 



on 



5A 



4tt 



-(F u " + 4tt K^l 



(50) 



(51) 



does not contribute to the variational principle J d 4 xSC = 0. As a result of the reduced variational principle, where 
the variations S and 8A U are arbitrary, we obtain the reduced Vlasov equation (|33[) and the reduced {macroscopic) 
Maxwell equations 



d 
dx^ 



4tt 



(52) 



from which we recover the reduced Maxwell equations (|42|) and (|43| . In addition, we find that the reduced charge 
conservation law (|4"Tj) follows immediately from the reduced Maxwell equations (f5"2")) . since the tensors F M1/ and K^" 
are antisymmetric. 



9 



2. Reduced Energy-Momentum Conservation Law. 



We now derive the energy-momentum conservation law from the reduced Noether equation SC = d^K^ associated 



with space-time translations x p 



dx^. The corresponding Eulerian variations {S,8A V ,8C) are 



S = PftSx* 1 
SA^ = 
5C = 



F M „ 5x 



5x v ) 



(53) 



After some cancellations introduced through the reduced Maxwell equations (|52p . we obtain the reduced energy- 

0, where the reduced stress-energy tensor is defined as 



momentum conservation law d^T^ 



47T 



^— F : F - (F^ + 4tt K^) F/ 

4 



jr--A u ) t 



(54) 



naturally includes the polarization and magnetization (j49|) . An explicit proof of the reduced energy- momentum 
conservation law is given in Ref. [24[. Note that, because of polarization and magnetization effects (^e ^ 0), the 
reduced stress-energy tensor (|54l) does not appear to be symmetric (while it may be so when evaluated explicitly). The 
physical meaning of this unsymmetrical Minkowski form [53] for the reduced stress-energy tensor has a rich history 
in physics (briefly discussed by Boghosian [TtJ and recently reviewed in Refs. [H, H^|) but its further discussion is 
beyond the scope of this paper. 



D. Push-forward Definitions of Polarization and Magnetization 

A complementary method for deriving the polarization and magnetization (|49p is provided by the push-forward 
method [24j ] . From the scalar- invariance relation l|27p , we construct the push-forward relation of fluid moments of an 
arbitrary function \ 011 particle phase space as follows. First, we transform the momentum average [x\ of an arbitrary 
phase-space function %: 

n[ X ] = Jd 3 pf X = J o?zT5\x-r) X 
= J d*zT T7 1 [S\x - r) x] 

= J d 6 zF S 3 (5i + p 6 - r)T 'fx, (55) 

where the extended Vlasov distribution (I34|) was used (n denotes the particle fluid density) and energy-time integra- 
tions were carried out. The reduced displacement (|32[) has components that are dependent on the fast space-time 
scales and components that are independent. Next, by expanding the right side of Eq. (|55p in powers of p e and 
integrating by parts, we obtain the push- forward relation 

n[ X ] ee n [T7 1 X ] A - V- (n [p £ T £ ~ 1 x ] A - V- \ [ PePe T £ ~ 1 X ] A + ■••), (56) 

where [• • -] A denotes an average with respect to the reduced Vlasov distribution F in reduced phase space (n denotes 
the reduced fluid density). In Eq. ([55]) . the terms in the divergence on the right side include the dipolar contribution 
(linear in p e ) and the quadrupolar contribution (quadratic in p t ). We now apply the push-forward relation (|56[) to 
derive dipolar expressions for the polarization and magnetization (where the quadrupolar and higher-order multipole 
contributions are ignored). 

First, we consider the case x — e (= T^x) so that the push-forward relation ([5^)1 yields the expression for the 
charge density 

p = p - v - (E e "t^] A + ■■■)' (57) 

where p = ^2 en defines the reduced charge density and higher-order multipole terms (e.g., quadrupole term that is 
quadratic in p e ) are not displayed. By comparing Eq. (|39| with the relation (|57j) between the particle charge density 
p and the reduced charge density p, we find the expression for the polarization (dipole-moment density) 

P ^ ^T(n[7r e ] A ) = ]T Jd 4 pn e T, (58) 
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where the reduced electric-dipole moment 



7T f 



e p e 



(59) 



is associated with the charge displacement induced by the phase-space transformation (f2l))) . Note that the fast-time 
average (tt 6 } = e (pj) of the reduced electric-dipole moment (|59[) may not vanish and thus may provide important 
polarization effects 0, H3] ■ 

Next, we consider the case x — ev = edx/dt, where the Lagrangian representation for the particle's velocity 
v = dx/dt is used. To obtain an expression for the push-forward T~ 1 (edx/dt), we use the reduced Vlasov operator 
d e /dt, defined in Eq. (|33|) . to obtain 



:T 



-i— 1 dx 

v = eT, — = e 
£ dt 



d e i 



dt 



dt 



(60) 



where we used Eq. (|31[) to obtain the last expression. The reduced (e.g., guiding-center) velocity d £ x/dt is independent 
of the fast space-time scales, while the reduced displacement velocity d e p c /dt exhibits both fast and slow space-time 
scales. Note that the fast-time-average particle polarization velocity d e (p e )/dt may not vanish [30j and thus may 
represent additional reduced dynamical effects (e.g., the standard polarization drift in guiding-center theory |6J) not 
included in the reduced velocity d e x/dt. 

We now replace the reduced moment [d e p e /dt] A in the first term of Eq. (|56[) by taking the time derivative of the 
polarization l|58p and using the reduced Vlasov equation ([55]) to obtain [2J] 



dP 



dp t 



dT 



E 



dePe 



dt 



dt 

A 



dt 



E 



e n 



d e x 
~~dt 



(61) 



where integration by parts of the reduced Vlasov equation 

d r x dT d £ p 



dt 



—r — = — 



dt 



dp dt 



dt 



d_ 
dp 



d e P 

dt 



was carried out to obtain the last expression. Hence, the push-forward relation (|56[) for x — edx/dt yields the 
expression for the current density [H| 



<9P 

~dt 

dP 

~dt 



c V x M 



e n 



Pe X 



d e x 
~dt 



vx E 



P e x 



dt 



(62) 



where J = ^ en [d e x/dt] A denotes the reduced current density, the second term represents the polarization current 
Jpoi = dP/dt, and the third term represents the magnetization current J mag = cV X M. By comparing Eq. (140]) 
with the relation (|62p , we obtain the magnetization 



1 d £ p e 



#>-E 



2 dt dt 

which is represented as the sum (for each particle species) of an intrinsic magnetic-dipole contribution 

deP e 



P, 



e 
27 



Pe X 



dt 



(63) 



(64) 



and a moving electric-dipole [23] contribution (n e X d e x/dt). We note that, since the relation n' = 7 n between the 
laboratory density n' and the rest-frame density n involves the relativistic factor 7 = 1/ ^/l — |v] 2 /c 2 , the standard 
moving-dipole magnetization contribution (3lj | to the polarization appearing in the charge density (|57j) comes from 



the last term in the relativistic c 



correction 



e 
2? 




e 




2^2 





d e x 


2 

+ 


d e /j e 


1 


dt 




dt 





e d e x 
? ~dt 



deP e 

dt 



(65) 
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where the first term yields a small relativistic correction to the electric dipole moment. We return to this relativistic 
magnetization effect in Sec. HVEl [see Eq. (|129[) ]. 

A simple example of the magnetization (|64|) is provided by guiding-center dynamics 0, [f| , where p denotes the 
gyroangle-dependent gyroradius vector (with d e p e /dt = Qdp/d6 to lowest order, where VL denotes the gyrofrequency) 
and the guiding-center intrinsic magnetization p gc = — /i b is defined in terms of the magnetic moment /i = m|v^ \ 2 /2B 
(an adiabatic invariant for guiding-center motion) and b = B/B. The guiding-center electric dipole moment on 
the other hand, is expressed as 7r gc = e(b/Q) Xv gc , where v gc denotes the guiding-center drift velocity (excluding 
the E x B velocity). 



E. Polarization and Magnetization 



We now summarize the main results of this Section. The polarization and magnetization 



P 

M 



E 
E 



<9E 
dB 



E " t 71 " 6 

E 



7T e d e X 

M e + — x — 

c dt 



can either be derived by variational method from the reduced potential 'J'e or by push-forward method from the 
reduced displacement p e = T^ 1 x — x. By combining these expressions, we obtain the relations between the reduced 
dipole moments (Tz e ,p e ) and the partial derivatives (d^ t /dE, d^ t /dB) of the reduced Hamiltonian: 



- <9* £ /<9E = 7r e 

- d^JdB = p e + (7r e /c) X dS/dt 
These relations emphasize the complementarity of the reduced variational and push- forward methods. 



(66) 



IV. OSCILLATION-CENTER DYNAMICS IN WEAKLY-MAGNETIZED PLASMAS 



The problem of the low-frequency oscillation-center Hamiltonian dynamics of charged particles in high-frequency, 
short-wavelength electromagnetic waves is the paradigm for applications of Lie-transform perturbation methods in 
plasma physics [l(| HH, EH El ■ Here, the fast wave space-time scales are removed asymptotically from the Hamiltonian 
dynamics by a time-dependent phase-space noncanonical transformation z a = (ct,x;w/c, p) — > z a = (c£,x; W/c, p), 
where x denotes the oscillation-center position, and (w, p = m v) denote the oscillation-center's kinetic energy- 
momentum. Note that, although the phase-space transformation is time-dependent, the time coordinate t remains 
unchanged by the transformation (i.e., the particle and oscillation-center times are identical). 

First, we decompose the electromagnetic potentials ($, A) in powers of the wave amplitude (represented by the 
ordering parameter e): 

<]> N ^^f;»i (67) 



' 71=0 V 

where the lowest-order background plasma is represented by the zeroth-order fields (<&o, Ao), while the primary wave 
fields (Ei,Bi) are represented by the first-order potentials (<&i,Ai) and the second-order (back-reaction) electro- 
magnetic fields (E 2 ,B 2 ) are represented by the second-order potentials (<J> 2 ,A 2 ). These second-order fields, which 
represents the plasma response to the external first-order fields, contain contributions that are weakly space-time de- 
pendent compared to the rapid wave space-time scales, and contributions with rapid second-harmonic wave space-time 
scales. 

The extended-phase-space Hamiltonian dynamics of charged particles in the electromagnetic fields (|67|) is expressed 
in terms of the extended Hamiltonian 

U = |p| 2 - w = Ho, (68) 
2 m 
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and the extended-phase-space Lagrangian 

T= p+ - (A + eAi + e 2 A 2 H ) • dx - [w + e ($ + e $1 + e 2 $2 H )] df 

= T + eT 1 + e 2 T 2 + (69) 

Here, the unperturbed Poisson bracket { , }o is obtained from To as 

( dF dG dF 8G\ e dF dG 

A. Eikonal Representation 

We now introduce the eikonal representation for the electromagnetic fields (f6"T|) . The previous work by Cary 
and Kaufman [TTI ] considered the case of an unmagnetized background plasma (with $ = = A ) and the case 
of a strongly magnetized plasma (with <!>o = 0). Here, we consider the intermediate case of a weakly-magnetized 
background plasma described as follows. 

First, we assume that the zeroth-order (background) fields = ($0, Ao) are weakly space-time dependent: 

Af 0) = Afafax), (71) 

where eo C 1 denotes the eikonal ordering parameter (warning: eo is not the permitivity of free space), so that the 
background electric and magnetic fields 

F(o)/iK = £ o F( ) M y(eoa;) (72) 

are ordered small; the weak-background-field (WBF) ordering (|T2"|) is consistent with a weakly-magnetized quasi- 
static background plasma. In particular, according to this ordering, the wave frequency u is considered much larger 
than the gyrofrequency f^o = eBo jmc of a charged particle moving in this weak magnetic field while the wavelength 
A = 2ir I It I 1 is considered much smaller than the particle's gyroradius po- Because of the mass dependence associated 
with this ordering, the weak-background-field (WBF) ordering (fT2"|) is especially appropriate for ions (because of their 
large masses), while the ordering might be invalid for electrons (because of their small masses). 
Next, the first-order wave fields are decomposed in terms of the eikonal representation: 

M = ( *O e «V" + ( ?* V* e/eo , (73) 



Ai J \A 1 J \Al 

where the wave eikonal-amplitudes (denoted by a tilde) are weakly-varying space-time functions and derivatives of 
the eikonal phase O(eor, eot): 



99 

eo 1 ^^) = fc^(e x) = (-w/c, k) (74) 



define the weakly- varying wavevector k and wave frequency to. Note that the first-order electromagnetic field has the 
eikonal-amplitude 

F(i) /JI / = i (k^ A {l)v - k v A(i)!j\ ■ (75) 
Lastly, the eikonal representation of the second-order fields is 

!O s (li) + (t:) e2,e/ " + (l) e " 2,e/ " ™ 

where the eikonal-averaged weakly- varying back-reaction fields ($2, A2) are generated by the second-order pondero- 
motive effects (as discussed below) while second-harmonic terms ($2, A2) are, henceforth, ignored (i.e., an eikonal- 
averaged expression involving second-harmonic terms appears at the fourth-order in wave amplitude, which falls 
outside the scope of our work). 
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B. Oscillation-center Transformation and Ponderomotive Hamiltonian 

As a result of the perturbation expansion (|6T[) in the electromagnetic potentials (<&, A), we see that the Hamiltonian 
exhibits explicit dependence on the eikonal phase and, thus, the particle dynamics exhibits fast space-time 
scales in addition to the slow space-time scales associated with the background plasma. 

Lie-transform Hamiltonian perturbation theory Q is concerned with the derivation of a new extended oscillation- 
center Hamiltonian 

H = T~ l H = (H Q + eHi + e 2 H 2 + •••)- w, (77) 
where the first-order and second-order oscillation-center Hamiltonians are 

Hi = Hi - d • dHo = - Gi • dHo, (78) 
H 2 = H 2 - G 2 -dH - - Gi -d(fli + Hi) = -G 2 -dH - - G x • dH x , (79) 

The last expressions in Eqs. (|78ll - (f79|) use the fact that the Hamiltonian perturbation terms H n (n > 1) are identically 
zero in the Hamiltonian (|68|> . The generating vector fields (Gi, G 2 , ■ ■ ■), on the other hand, are determined from the 
requirement that the oscillation-center phase-space Lagrangian 

r = T^r + ds = T , (80) 

retains its unperturbed form. From this requirement, we obtain 

GI = {Si, z a } + - c A {l) ^ z a } , (81) 
G a 2 = {S 2 , z a } + tA mit {x», z a } - ± {Si, x»} F (1)/1V {x», z a } 0l (82) 

where the scalar fields (Si, S 2 , ■ ■ ■) are determined from Eqs. (|T8|) - (fT9")) by requiring that the oscillation-center Hamil- 
tonian (|77|) be eikonal-phase independent. 

1. First-order analysis. 

From Eqs. ijTSjl and (f8"Tj) . we obtain the first-order Hamiltonian 

Hi = - {Si, H Q } - - c A {x)ll H }o = - ^ + e($i - ^-Ai) = 0, (83) 

where the first-order oscillation-center Hamiltonian Hi must vanish since (<£>i, Ai,Si) are all explicitly eikonal-phase 
dependent. By inserting the eikonal representation (|73|) into Eq. (|83|) . we obtain the eikonal equation for Si: 

-iu'Si + e e (Eo + ^xBo) • = e ($i - ^-Ai) (84) 

where us' — lu — k • v is the Doppler-shifted wave frequency, and the WBF eikonal ordering (|72)) has been used. Next, 
by substituting the WBF-eikonal expansion 

Si = Sio + eoSu + • • ■ (85) 
into Eq. we obtain the lowest-order eikonal amplitude 

Sio = -,(®i - --Ai), (86) 

and the first-order WBF-eikonal correction 

Sn = - — (e + -xB ) -l = -<F -4» (87) 
oj' \ c ) uj' 
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where the first-order spatial displacement £ is defined in terms of the lowest-order eikonal amplitude [1J 

~: dS 10 e 



dp 



TflLO 



Ei + - x Bi 

c 



(88) 



Note that the first-order spatial displacement £ also satisfies the eikonal relation imio' £ = ikSio +eAi/c, which can 
be expressed as 



d $, a 
- m —j- = V610 + - Ai, 
dt c 

where do/dt = d/dt + v • V and u>' = — £q 1 doQ/dt. 

The phase-space transformation (|26[) can thus be expressed explicitly as 

x = x — edS w /dp + ■ ■ ■ = x — e£ + ••• 

p = p + e (VSio + eAi/c) H = p - emd $,/dt + 

w = w — e (dSw/dt — e<I>i) + • • • = w — emv ■ do^/dt 
Hence, to first order in e, the oscillation-center momentum coordinate is 

doC 



m v — e 



dt 



and the oscillation-center energy coordinate is 



dt 



~dt' 



2m 



(89) 



(90) 



which satisfies the unperturbed constraint 7io = |p| 2 /2m — w = for the oscillation-center kinetic energy- momentum 
coordinates. 



2. Second-order analysis. 



Next, from Eqs. ([79]) . {82]) , and ([83]) . the expression for the second-order Hamiltonian is 



ff 2 



e $o A- 



e /„ v „ 
- E1 + -XB1 

z V c 



agi 

<9p 



(91) 



where the right side has both eikonal-dependent and eikonal- independent terms. To derive the eikonal- independent 
terms that define the second-order oscillation-center Hamiltonian H2, we perform an eikonal-phase average of the 
right side of Eq. (|9Tj) to obtain 



H 2 = e $0 - 



(( E -7-.)-f). 



(92) 



where the contributions of the second-harmonic eikonal-amplitudes ($2, A2,S , 2 = 5*2) vanish. 
Using the definition ([88]) , the lowest-order contribution from the second term in Eq. ([92]) is 



Ei 



x Bi 



r ' dp 



(93) 



which determines the standard ponderomotive Hamiltonian [Tl(. Using the definition ([57]) . on the other hand, the 
first-order WBF correction is 



(( 



E 1 + ^xBi).fH 
c J op 



dp 



(e + ^ X B ) -7f 2 - B -7I 2 , 



dp 



• »' F n X 



(?xf) 



(94) 
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where we introduced the definitions for the ponderomotive electric dipole moment 

tF 2 = ekx (ilxT) , (95) 
and the ponderomotive magnetic dipole moment 

p 2 = - J (i£xT) ■ (96) 



Note that the intrinsic ponderomotive magnetic-dipole moment p 2 is accompanied in Eq. (|94p by the ponderomotive 
moving-electric-dipole moment 7r 2 X v/c 27]. As indicated above, however, the moving-magnetic-dipole contribution 
(v/c) X p 2 to the ponderomotive polarization appears as a relativistic correction, which is ignored in the present 
nonrelativistic treatment. 

3. Ponderomotive Hamiltonian. 

By combining Eqs. (193)) - ([94)1 . the second-order oscillation-center Hamiltonian is defined by the eikonal-averaged 
terms on the right side of Eq. (|9"Tj) : 



H 2 = e ( $ 2 - - • A 2 ) + muj' 2 |£| 2 - e 



E • 7T 2 + B • I £t 2 + 7T 2 X - 



c 



c 



e[ $ 2 - --A 2 ) + * 2 , (97) 



where the reduced potential ^2 = ^20 + eo ^21 contains the ponderomotive potential ^ 2 o and its first-order (WBF) 
correction ^21, which includes the ponderomotive polarization and magnetization effects (|95|) -(|96 |) . 

Lastly, by using the ponderomotive dipole moments ([M]) and we can immediately write down expressions for 
the ponderomotive polarization and magnetization 

(F„ B.) s E / <e? F ft + W 2 X = - £ / ^ !• ^) , (98, 

which satisfy the relations (|66p between the variational and push- forward methods. Note that the polarization and 
magnetization 

P = e Pi + e 2 P 2 + ••• and M = e M a + e 2 M 2 + • ■ • (99) 

have first-order contributions that vanish upon eikonal-phase averaging (i.e., (Pi) = = (Mi)), while the eikonal- 
averaged contributions are second-order in wave amplitude: (P) = e 2 P2 and (M) = e 2 M2. 

C. Oscillation-center Pull-back Operator 

The relation between the particle Vlasov distribution / and the oscillation-center Vlasov distribution F is expressed 
in terms of the pull-back operator T e : 



G 2 ■ d,F + - Gi ■ d (Gi ■ d,F) 



/ = T e F = F + e Gi • dF + e 2 

= fo + e h + ( 2 h + ■ ■ ■ , (100) 
where the vectors (Gi, G2, • ■ ■) are defined in Eqs. (f8Tj) and (|82|) . 



1. First-order pull-back operator. 

Since the oscillation-center Vlasov distribution F is independent of the fast wave space-time scales (i.e., it is 
independent of the eikonal phase 0), the eikonal-phase averaged part (/1) of the first-order particle Vlasov distribution 
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fi = Gi • dF vanishes: (/i) = (since (Gi) = 0). Its eikonal-phase dependent part /i = /io + eo /n + • • •, on the other 
hand, is expressed in terms of the lowest-order contribution 



r / t dF ie r~ v 

f 1Q = imu'£- — = -— (E1 + -XB1 

Op LU' V C 



dF 
dp ' 



F (ik-a d 3 p, 



so that the zeroth velocity-moment is 

j fio d 3 P = -ij F A . d 3 p = - 

from which we obtain the first-order polarization eikonal-amplitude 

?i = E e / ^^ 3 p- 

The first velocity- moment, on the other hand, is 

vf w d 3 P = -i J F A. (w£v) d 3 p = - J F [(ik-tj v + tw'£ 



d 3 p 



F 



-i u; £ + ik X ( £ X v 



(101) 



(102) 



(103) 



(104) 



where the first term in the last expression is associated with the first-order polarization current, while the second 
term is associated with the first-order magnetization eikonal-amplitude 



Mi = e / F t x ~ d3 p 



(105) 



Note that only the moving electric-dipole contribution appears here since the intrinsic magnetization must be quadratic 
inf. 



2. Second-order pull-back operator. 



The second-order particle Vlasov distribution / 2 = G2 • dF + i Gi • d(Gi • dF) has an eikonal-phase independent 
part expressed as 



</20> - ^ 



~ (A 2 + Re(Bi X £*)) F + mV 2 Re (| £*) • A 



(106) 



where only the lowest-order WBF terms are kept. Note that this second-order particle contribution has a vanishing 
zeroth velocity-moment 

(107) 



</ 20 ) d 3 p = 0, 

since Eq. (|106|) is expressed as a momentum-space divergence. However, the first velocity-moment 



v (/20) d 3 p = - d 3 p 



— ( eA 2 



mc mc 



— Re (Bi x £ 

mc V 



— * dF 



F 



— A 2 + L0'k- ££ +£ £ 

TOC 



(108) 



is non-vanishing, where we used the following identities in obtaining the last expression: d p (mLu' £) = £k and 



Re 



(Bi x£*) 



moj 



|£| 2 k - k-Re(£ I) 



which follows from the relation (e/c) Bi = — too/ k X £ obtained from Eq. Using the lowest-order second-order 

oscillation-center Hamiltonian (|QT[) . we obtain the identity 



v </ 20 ) d 3 p = / F A 2 

mc 



-r- 9* 



(109) 
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Lastly, the second-order kinetic energy 

/ d3p E) {ho) = / T (* 2 ° ■ ^ 2 '^) ^ (no) 

is naturally expressed in terms of the oscillation-center distribution F and the (lowest-order) second-order oscillation- 
center Hamiltonian ([W)) . 



D. Ponderomotive Variational Principle 

The variational formulation for the exact Vlasov-Maxwell equations was presented in Sec. HH The covariance of 
the Vlasov part of the action functional has been used in Ref. [32| to derive a variational formulation of the nonlin- 
ear gyrokinetic Vlasov-Maxwell equations, which describe the turbulent evolution of low-frequency electromagnetic 
fluctuations in a magnetized plasma with arbitrary geometry . 



1. Linear polarization and magnetization. 



The eikonal form of the first-order Maxwell equations can be expressed in terms of the eikonal-averaged second-order 
Lagrangian density 



1 

-in 
1 



Eil 



Bil 2 ) - E / d "P T 



e $ 



■•A; 



20 



-Re(Ex.Dj - Bx-Hl) - £ 



d p F ( $2 - - • A a 



(111) 



where the last expression is obtained by writing the ponderomotive potential i n terms of the first-order polarization 
and magnetization (|103|) and (|105[) as 



J d 3 pF * 20 = - Re (PI -Ei + • Bi 



which displays the standard relation between the ponderomotive Hamiltonian and the oscillation-center polarization 
and magnetization [H, [34[ . 

We define the first-order eikonal-amplitudes for the macroscopic electromagnetic fields [13] 



Hi 



4tt -J^- = Ei + 4ttPi, 



47T 



dCw 
dB\ 



Bi - 4ttMi, 



(112) 
(113) 



in terms of the first-order eikonal-amplitudes for the polarization and magnetization (|103| and (|105|) . Hence, by using 
these definitions, the first-order Maxwell equations become 











dC 20 _ dE| _9£2o = . k Di 
d$l ' d$l 9E* 4tt' 
dC 20 _ 8EI dC 20 , dBl 3C 2 o 



= _.w Di _ jkx Hi 
dAl dAl d~El dAl dB{ c 4tt 4tt ' 

where the eikonal amplitudes for the first-order charge-current densities are 

-iuj \ + ikx X v 



(114) 
(115) 



Pi 
Ji 



E 



d 3 p fw 



E 



d 3 pF 



-ik-Pi 

-iwPi + ike X Mi 



(116) 



which are consistent with Eqs. (|103p and (|105p . Note that Eq. (| 1 1 5|) is consistent with Eq. (|114[) . i.e., the dot product 
of Eq. (|115|) with the wave vector k yields Eq. (|114|) . 
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2. Nonlinear polarization and magnetization. 

The derivation of the lowest-order expressions for the eikonal-averaged second-order charge and current densities, 
and the eikonal corrections corresponding to the ponderomotive polarization and magnetization, proceeds with the 
introduction of the second-order eikonal-averaged Vlasov-Maxwell action functional 



47T 



Re (Ei • T>1 - Bi • Hj) + e (E • D 2 B • H 2 
e fd 3 pF($ 2 - --A 2 j. (117) 



We first note that, since the second-order eikonal-averaged Vlasov Lagrangian density (|117| is independent of the 
background scalar potential 4>o, we find 

BC 

**> s - oi = °' ^ 

and, thus, the second-order oscillation-center eikonal-averaged charge density must be generated by ponderomotive 
polarization effects [see Eqs. (|121|) - (|122|l ] to first order (in eq) in eikonal analysis. This result is consistent with the 
fact that the zeroth-moment (|107j) of the second-order eikonal-averaged particle Vlasov distribution (/20) vanishes. 

The second-order ponderomotive (eikonal-averaged) polarization and magnetization are obtained from the action 
functional (|117[) as 



<9E 4?r 



«■ - -° 'ft + 1 - E / + ' 2 x ') ^ (120) 

where W2 and 7*2 are defined in Eqs. (|95p and (|96[) . The second-order ponderomotive charge and current densities are 

p 21 = -V-P a , (121) 
J21 = y + cVxM 2 , (122) 

respectively. 

Lastly, the second-order eikonal-averaged Maxwell equations can also be written as 

V-D 2 = and V X H 2 - = — J 20 , (123) 

c at c 

which requires J20 to be divergenceless [TTj] , as is also required for the eikonal-averaged second-order charge and current 
densities 7j 2 = P20 + e o P21 an£ l J2 = J20 + £o J21 to satisfy the eikonal-averaged second-order charge conservation law 

^ + V-J 2 = 0. (124) 

Hence, the condition (|118[) immediately implies that V-J20 = [ll[, while the second-order ponderomotive charge 
and current densities (|12ip - (|122|) explicitly satisfy Eq. (|124p . 

We have, thus, achieved the purpose of this Section, which was to show how the ponderomotive polarization and 
magnetization ([55)) can be self-consistently derived from a variational principle (|117[) for the oscillation-center Vlasov- 
Maxwell equations. We note that ponderomotive polarization and magnetization effects in magnetized plasmas have 
also been investigated in Refs. [3, HH by using a Low-Lagrangian-type variational formulation. 



E. Ponderomotive Push-forward Derivation 



We now proceed with the complementary derivation of the ponderomotive polarization and magnetization by the 
ponderomotive push-forward method. By inserting the oscillation-center transformation generated by Eqs. {5TJ and 
(|82f into the reduced displacement (|32|) . we obtain the oscillation-center displacement 



p e ee e£ + - [€.VC + m§-|i) - e 2 G5 + (125) 



2 V dt °P 
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where Eq. (|89[) was used and the first-order WBF corrections are omitted in what follows. While the second-order 
term = — d p S2 appears on the right side of Eq. Q125p , it does not contribute to the relations derived below because 
its contributions are of higher order in e. 



1. Linear polarization and magnetization. 



The linear displacement p f = e£ is eikonal-phase dependent and its amplitude is p t = e£. Hence, the push-forward 
formula (158p yields the linear polarization 



Pi = > en £ = — > n — — — 
^ L J ^ 9E* 



(126) 



The push- forward formula (|62p . on the other hand, yields the linear expression for the eikonal-dependent current 
Ji 



-iw Pi + ik X 
so that the linear magnetization is 



E 



£ x v 



■iwPi + ike X Mi 



c 



E 



dBl 



(127) 



where we used the lowest-order expression x e = v for the reduced particle velocity. Note that the first-order magne- 
tization (|127|) only has a moving electric-dipole contribution (i.e., Jjl 1 = 0) since the intrinsic magnetization (|64|) is 
quadratic in p £ . 



2. Ponderomotive polarization and magnetization. 

Ponderomotive polarization and magnetization effects enter at second order through eikonal-phase averaging. We 
begin with the eikonal-averaged reduced displacement 

= e 2 kx (*txT) = e 2 (^^2), 

from which we recover the ponderomotive electric-dipole moment (|95|). Hence, the push- forward relation (|58|> yields 
the second-order ponderomotive polarization 



p 2 = e - 2 e = E = - E » ( 128 ) 

We note that the ponderomotive magnetization contribution to the ponderomotive polarization is obtained from the 
relativistic correction (|65|) as 

«< ft T7'(7-i)>- ?(*(*-^)) + - 

s I X (« 3 p 2 ) , (129) 

where d e x/dt = v (to lowest order in e) and we have omitted the relativistic correction (|v| 2 /2c 2 )(p e ) to the pon- 
deromotive polarization (|128p . Furthermore, we note that the quadrupolar contribution to the ponderomotive po- 
larization (|128p is of the form e 2 V • Q 2 , where the ponderomotive electric quadrupole moment (tensor) is defined as 

Q 2 = ^2 en Re([££ ]), which is one order higher than the ponderomotive polarization in the eikonal analysis (i.e., 
V • Q2I ~ eo |I*2 1) and ponderomotive quadrupolar polarization is therefore omitted. 

Next, we derive an expression for the second-order ponderomotive magnetization. First, the moving-dipole contribu- 
tion is ((e p e X x c )) = e (p e ) X x e = e 2 W 2 X v to lowest order in e. Second, the intrinsic ponderomotive magnetization 
is 

((!^t)) = 8 ^ 



e 2 Re 



ik- £ £ 



9£_ 
dp 
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where we used (p e , d e pj dt) — > (e£, — iea/£). The push-forward relation (I63|) yields the second-order ponderomotive 
magnetization 



M 2 = e- 2 



P e x o 



1 d e p £ 



2 dt 



d e x 
~dT 



E 



/X 2 + 7T 2 X — 

C 



E 



9B 



(130) 



Equations (|128[) and (I130|) represent the main results of this Section. They exhibit the complementarity of the 
push-forward and variational methods for the case of the oscillation-center dynamics of charged particles in a weakly- 
magnetized background plasma. 



V. SUMMARY 



From a historical perspective, the derivation of oscillation-center Hamiltonian dynamics provided an ideal battle- 
ground for the application of Lie-transform perturbation methods in plasma physics, in which the Berkeley School 
of Plasma Physics, led by Allan N. Kaufman and his collaborators, played a fundamental role. In the spirit of the 
Berkeley School, the present paper demonstrated the interconnectedness of the Lie-transform perturbation method 
with the variational approach. 

The variational formulation (|117|) of the oscillation- center Vlasov-Maxwell equations provides a simple description 
of self-consistent ponderomotive polarization and magnetization effects in weakly-magnetized plasmas, derived by vari- 
ational derivatives with respect to the background electromagnetic fields (Eo, Bo). Although the treatment presented 
here does not include relativistic effects, the variational formalism is naturally suitable for relativistic generalization. 

Ponderomotive effects for unmagnetized plasmas [ill ] can be easily recovered by setting Ao = in the final expres- 
sions. Ponderomotive effects in strongly-magnetized plasma, on the other hand, require a two-step transformation: 
first, a transformation from particle to g uiding-center phase space is introduced [J] to be followed by a transformation 
to oscillation-center phase space fToL I 111 ]. Lastly, low-frequency ponderomotive polarization and magnetization effects 
in magnetized plasmas arc described in terms of the nonlinear gyrokinetic Vlasov-Maxwell equations (e.g., see Ref. [8j) 
or the nonlinear electromagnetic drift-fluid equations [20j. 
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